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Abstract 

We study collective excitations in a helical electron liquid on a surface of three-dimensional topological insulator. 
Electron in helical liquid obeys Dirac-like equation for massless particles and direction of its spin is strictly 
determined by its momentum. Due to this spin-momentum locking, collective excitations in the system manifest 
themselves as coupled charge- and spin-density waves. We develop quantum field-theoretical description of spin- 
plasmons in helical liquid and study their properties and internal structure. Value of spin polarization arising in the 
system with excited spin-plasmons is calculated. We also consider the scattering of spin-plasmons on magnetic 
and nonmagnetic impurities and external potentials, and show that the scattering occurs mainly into two side 
lobes. Analogies with Dirac electron gas in graphene are discussed. 
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1 Introduction 

Topological insulator is a new class of solids with non- 
trivial topology, intrinsic to its band structure. Theoreti- 
cal and experimental studies of topological insulators 
grow very rapidly in recent years (see [1,2] and refer- 
ences therein). Three-dimensional topological insulators 
are insulating in the bulk, but have gapless surface states 
with numerous unusual properties. These states are 
topologically protected against time-reversal invariant 
disorder. When gap is opened in surface states by time- 
reversal or gauge symmetry breaking, a spectacular mag- 
netoelectric effect arises [3,4]. 

Recently a "new generation" of 3D topological insula- 
tors (the binary compounds Bi 2 Se 3 , Bi 2 Te 3 and other 
materials), retaining topologically protected behavior at 
room temperatures, were predicted and studied experi- 
mentally [5-7]. Band structure of the surface states of 
these materials contains a single Dirac cone, where elec- 
trons obey 2D Dirac equation for massless particles. 
Direction of electron spin in these states is strictly 
determined by their momentum, so these states can be 
called as "helical" ones. Surface of topological insulator 
can be chemically doped, forming charged helical liquid. 
The spin-momentum locking leads to interesting 
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transport phenomena including coupled diffusion of 
spin and charge [8], inverse galvano-magnetic effect 
(generation of spin polarization by electric current) [9] 
and giant spin rotation on an interface between normal 
metal and topological insulator [10]. The spin-momen- 
tum locking offers numerous opportunities for various 
spintronic applications. 

Collective excitations (plasmons) in helical liquid on 
the surface of topological insulator was considered in 
[11]. It was shown that due to spin-momentum locking 
responses of charge and spin densities to external elec- 
tromagnetic field are coupled to each other. Therefore 
the plasmons in the system should manifest themselves 
as coupled charge- and spin-density waves and can be 
called "spin-plasmons". In [12] application of spin-plas- 
mons in spin accumulator device was proposed. Also 
surface plasmon-polaritons under conditions of topolo- 
gical magnetoelectric effect were considered in [13]. 

Properties of the states on a surface of 3D topological 
insulator are similar to those of electrons in graphene. 
Graphene is unique 2D carbon material with extraordinary 
electronic properties [14-16]. Its band structure contains 
two Dirac cones with electrons behaving as massless Dirac 
particles in their vicinities. Graphene is a perspective 
material for nanoelectronics due to large carrier mobilities 
at room temperature. Electronic interactions and collective 
excitations in graphene have been extensively studied (see 
[17] and references therein). In particular, the properties 
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of plasmons [18,19] and hybrid plasmon-photon [20] and 
plasmon-phonon [21,22] modes were investigated theore- 
tically and experimentally. It was realized recently that gra- 
phene is a fertile ground for quantum plasmonics [23] due 
to very small plasmon damping. 

In this article, we develop quantum field-theoretical 
formalism to describe plasmons in graphene and spin- 
plasmons on a surface of 3D topological insulator based 
on random phase approximation (RPA). Problems of 
excitation, manipulation, scattering and detection of sin- 
gle plasmons can be conveniently considered using this 
approach. Thus, this approach can be especially useful 
for problems of plasmon quantum optics and quantum 
plasmonics. We use our approach here to study internal 
structure and properties of spin-plasmons in a helical 
liquid. 

The rest of this article is organized as follows. In Sec- 
tion 2, we present a brief description of electronic states 
on a surface of topological insulator and in graphene. 
Next we develop the original quantum field-theoretical 
description of plasmons in Dirac electron gas in Section 
3 and apply it further to study their properties. We con- 
sider internal structure of plasmons in Section 4 and 
important consequences of spin-momentum locking in 
Section 5. Scattering of plasmons on impurities and 
external potentials is considered in Section 6, and Sec- 
tion 7 is devoted to conclusions. 

2 Dirac electrons 

The low-energy effective Hamiltonians for electrons in 
helical liquid [24] and graphene [16] are 

helical liquid: H 0 = v^[p x a y - p y a x ), (1) 

graphene: H 0 = v ¥ {p x cr x + p y cr y ), (2) 

where the Fermi velocities of electrons v F are 6.2 x 
10 5 m/s for the topological insulator Bi 2 Se 3 and 10 6 m/s 
for graphene; the Pauli matrices o x and o y act in the 
spaces of electron spin projections (helical liquid) or 
sublattices (graphene). The eigenfunctions of (l)-(2) can 
be written as e lv ' x \f py )l\fS, where S is the system area 
and [f p7 ) is the spinor part of the eigenfunction, corre- 
sponding to electron with momentum p from conduc- 
tion (y = +1) or valence (y = -1) band: 

helical liquid: \f PY ) = — .^ p/2 J , (3) 

graphene: \f py ) = — ^ ^ p/2 j , (4) 

where 0 p is a polar angle of the vector p (here and 
below we assume h = 1). Another difference between 



helical liquid on a surface of topological insulator and 
electron liquid in graphene is additional fourfold degen- 
eracy g = 4 of electronic states in graphene by two spin 
projections and two nonequivalent valleys. 

The value of electron spin in helical liquid and of 
pseudospin in graphene in the state \f p7 ) is 

helical liquid: (f py \cr\f py ) = y[zx p], (5) 

graphene: (f py \a\f py ) = yp, (6) 

where p and z are unit vectors directed along the 
momentum p and the z-axis, respectively. We see that, 
in helical liquid, the spin of electron lies in the system 
plane and makes an angle 90° (in counterclockwise 
direction in the conduction band and inversely in the 
valence band) with its momentum. In graphene, the 
sublattice pseudospin of electron is directed along its 
momentum in conduction band and opposite to it in 
the valence band. Physically, a definite direction of the 
pseudospin in the system plane corresponds to definite 
phase shift between electron wave functions on different 
sublattices. 

A starting point for quantum field-theoretical consid- 
eration of plasmons on the surface of topological insula- 
tor and in graphene is the many-body Hamiltonian of 
electrons with Coulomb interaction between them: 

H =J2 %PY a+ PY a PY + ^HI212 v <i 

PY qpp' YiYi y[Yi 

(7) 

x (fp+q,K{ 1/pKi ) (/p'-q,K2 \fp'Yi ) 

x a p+q,W a p'-q^2 a ^' yi dpYl ' 

where IS the destruction operator for electron with 
momentum p from the band y, £ p7 = yv^p-^ is its kinetic 
energy measured from the chemical potential ^ and V q 
= 2ne 2 lsq is the 2D fourier transform of Coulomb inter- 
action potential screened by surrounding 3D medium 
with a dielectric permittivity e. 

3 Description of plasmons 

To investigate the properties of plasmons in Dirac elec- 
tron gas, we develop the equation-of-motion approach 
similarly to the original works on plasmons in usual 
electron gas [25]. We treat a plasmon as a composite 
Bose quasiparticle, consisting of electron-hole pairs with 
the common total momentum q. Thus, the creation 
operator for plasmon with momentum q can be written 
as: 

Qq = a l+qy a py (8) 

PYY' 
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Here the coefficients Cp q K are the weights of intraband 
(y = Y) and interband (y ' * y) single-particle transitions, 
contributing to the wave function of plasmon. 

The plasmon creation operator should obey Heisen- 
berg equation of motion 



[h,q:i 



(9) 



where Cl q is plasmon frequency. We start from equa- 
tion of motion for single electron-hole pair, which can 
be derived using (7): 

[H,al +qy a VY \ 



y" 

(fpy \fp+q',y" ) a p+q,y' a p+q',y" ) 



(10) 



where we have introduced the Fourier transform of 
the density operator: 



Pq ~ tfp+q,Y'\fpy) a p+q l y' a py - 

pyy' 



(id 



The right-hand side of the Equation (10) contains pro- 
ducts of four fermionic operators. To reduce it to pro- 
ducts of two fermionic operators we use the RPA; its 
applicability will be discussed below. According to RPA 
[25], the operator products in the last two lines of (10) 
can be replaced by their average values in ground state | 
0) of the system: 



a p'y ,a v>y 



(0 ApyflpJ 



0) 



: $pp'$yy' i 



l py 



(12) 



where n p7 is the occupation number for electronic 
states with momentum p from the band y For electron- 
doped Dirac liquid at T = 0 (see also the remark at the 
end of this section), we have n p+ = ®(j?F - |p|)> = 1> 
where p ¥ - ^/v F is the Fermi momentum. In the case of 
hole doping, all characteristics of plasmons are the same 
due to electron-hole symmetry. 

Thus, the equation of motion (10) for electron-hole 
pair with taking into account the RPA assumption (12) 
takes the form: 



[H f a p+q, y > a py ] ~ (£p+q,x' £py ) a p+q,y' a py 



"p+q,K 

+ ~7T Pqifpylfp+qy) i n py ~ n p+q,K')- 



(13) 



Combining the definition of plasmon creation opera- 
tor (8) with the equation of motion for plasmon (9) and 
single electron-hole pair (13), we obtain the system of 
equations for the coefficients Cp q y : 



+ £py %p+q,y')Cpq - (fp+q,y'\fpy) 

x ^p' T qifp'r\fp'+q,T')i n p'T ~ n p'+q,x')- 
p'xx' 



(14) 



Introducing infinitesimal damping into the plasmon 
frequency and denoting 

^q = y I] C pq (/p*l/p+q^H>V " TWO/ (15) 

pxx' 

we can find the plasmon wave function from (14) as: 
{fp+q,y' \fpy)Nq 



Q 



pq 



+ Hpy ~ %p+q,y' + ^ 



(16) 



The normalization factor N q can be determined from 
the commutation relation for plasmon operators, which 
should be satisfied on the average in the ground state: 



<0|[Qq,Qq']|0) 

Substituting (8) in (17), we get 

T, D r'r = !' 

yy' 



(17) 



(18) 



where the total weights of intraband (D++) and inter- 
band (D + _ + D_ + - 1 - D ++ ) electron transitions, contri- 
buting to the plasmon wave function (16), are: 



n 



p+q,y' 



(19) 



To find the plasmon frequency Cl q , we can substitute 
(16) in (15), arriving to the standard RPA equation (see 
it for the helical liquid [11] and for graphene [18,19]): 



i-v«nfa, n«) = o, 



(20) 



where the polarization operator of Dirac electron gas 
is introduced: 



pyy' 



(21) 



lpy 



sp+q,y 



+ iS 



The degeneracy factor g is 1 for helical liquid on the 
surface of topological insulator and 4 for graphene. The 
angular factors (f p +q,f\fpy) are specific to helical Dirac 
electrons and arise in (21) as a result of summation over 
spinor components of electron wave function. 

The RPA becomes exact in the limit of small values of 
dimensionless parameter r s , defined as a ratio of charac- 
teristic Coulomb interaction energy to kinetic energy. 
For the gas of massless particles, r s does not depend on 
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electron density and equals to r s = e 2 /sv ¥ (effectively s = 
(si + s 2 )/2 when 2D electron layer is surrounded by two 
dielectric half-spaces with permittivities s x and e 2 ). For 
Bi 2 Se 3 , r s = 0.09 with e = 40 for dielectric half-space, and 
applicability of the RPA is well established (the value of r s 
for another material Bi 2 Te 3 is close to that for Bi 2 Se 3 ). In 
the case of graphene, r s is not small, but applicability of 
the RPA can be established due to smallness of the para- 
meter 1/g, leading to selection of bubble diagrams (see 
[16] and references therein, and also the work [26]). 

Maximal achievable amounts of doping of helical 
liquid on a surface of Bi 2 Se 3 are \a ~ 0.3 eV [24], there- 
fore at room temperature it can be degenerate electron 
liquid. Hence we assume that T = 0 in all calculations 
below. 

4 Wave function of plasmon 

Plasmon dispersions Q q , calculated numerically from (20)- 
(21) at various r s , are plotted in Figure 1. We also present 
results for suspended graphene with r s = 8.8 (s = 1) and 
for graphene embedded into Si0 2 environment with r s = 
2.2 (s = 4). For convenience, the degeneracy factor g = 4 is 
included here into r s , since the plasmon dispersion in RPA 
(20) depends only on the combination gr s . 

It is seen that for small values of r s (the case of topo- 
logical insulator) the plasmon dispersion law approaches 
very closely to the upper border of the continuum of 
intraband single-particle transitions (co <v ¥ q). On the 
contrary, at moderate and large r s , plasmon has well- 
defined square-root dispersion in the long-wavelength 



range. At q ~ p ¥ , the plasmon enters the continuum of 
interband single-particle transitions (co > 2(A - v ¥ q). 
Inside the single-particle continuum, energy and 
momentum conservation laws allow energy transfer 
between plasmon and single-particle excitations, so the 
plasmon acquires a finite lifetime. 

The internal structure of plasmons in Dirac electron 
gas can be characterized by total weights (19) of intra- 
and inter-band transitions in its wave function. The 
weight D ++ of intraband single-particle transitions is 
plotted in Figure 2. It is seen that the undamped plas- 
mon for all values of parameter r s consists mainly of 
intraband transitions, thus an influence of valence band 
on its properties is rather weak. When the plasmon 
enters the single-particle continuum, inter- and intra- 
band transitions start to contribute almost equally to its 
wave function, but the plasmon undergoes strong 
Landau damping. 

The detailed picture of the plasmon wave function in 
momentum space Cp q K > showing the distribution of con- 

I I 2 

tributions of intraband C + tl or interband 



pq 



\C + ~ 



Q 



pq 



electron-hole pairs into wave function of 

plasmon with given momentum q is presented in Figure 
3. The results are calculated for q = 0Ap F at r s = 0.09 
and 8.8. For other values of plasmon momentum q, the 
results are qualitatively similar. It is seen that for small 
values of r s (the case of topological insulator) the distri- 
bution of intraband contributions is very sharply peaked 
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Figure 1 Plasmon dispersions. Plasmon dispersions at different 
values of r s . Continuums of intraband {co <v f q) and interband {co > 
2jj - v f q) single-particle excitations are shaded; p F = \xN v is the 
Fermi momentum, \x is the chemical potential. 
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Figure 2 Weight of intraband transitions in plasmon wave 
function. Weight D ++ (19) of intraband electron-hole transitions in 
plasmon wave function at different values of r s . 
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in the forward direction, whereas the contribution of 
interband transitions is negligible. It can be understood 
from the reason that the plasmon dispersion is very 
close to the single-particle continuum (see Figure 1) and 
thus the plasmon itself behave almost as single intra- 
band electron-hole transition. 

At large values of parameter r s (the case of suspended 
graphene), the broad range of electron intraband transi- 
tions in momentum space contributes to plasmon. The 
weight of interband transitions is small but not negligi- 
ble. Contributions of interband transitions form two 
side lobes, since the angular factor (f p+qi ,±\f p +) suppresses 
interband forward scattering. 

5 Charge- and spin-density waves 

When a spin-plasmon is excited in the helical liquid, 
anisotropic distribution of electron-hole pairs of the 
type, depicted in Figure 3, arises. This distribution is 
shifted towards the plasmon momentum q. Due to the 
spin-momentum locking, the system should acquire a 
total nonzero spin polarization, perpendicular to q. A 
similar situation occurs in the current-carrying state of 
the helical liquid, which turns out to be spin-polarized 
[9]. 

Average of spin polarization in the state |lq) = Qq|0) 
with a single plasmon excited is (s) = (l q |o72|l q ) and 
can be calculated using (8) as: 



vyy' t (22) 

~~ Cpq (f pT |cr|/ py )j ^Cpq^ (%y — H p +q, y '). 

If q is parallel to e x , only the j-component of (s) is 
nonzero. Its dependence on q for undamped plasmons 
at various r s is plotted in Figure 4. At large enough r s , 
the spin polarization of the system is comparable with 
the whole spin of a single electron. 

Note that in the case of graphene the isospin polariza- 
tion of the system appears instead of spin polarization. 
Isospin polarization corresponds to nonzero average 
phase shift between wave function of electrons on differ- 
ent sublattices. On the contrary, in pseudospin-unpolar- 
ized state this shift is zero on the average. 

In ordinary electronic system without spin-momentum 
coupling, plasmon manifests itself as charge-density 
wave. In helical liquid on surface of topological insula- 
tor, due to the spin-momentum locking, spin-plasmon 
manifests itself as coupled charge- and spin-density 
waves. These waves, accompanying spin-plasmon with 
the momentum q, can be characterized by correspond- 
ing spatial harmonics of charge- (11) and spin-density 

S q = \ J2 tfv+W \ G \fx>Y) a l+q,y> a VY (23) 

PYY' 
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0.8 1 

Figure 4 Spin polarization of the system with an excited spin- 
plasmon. In-plane transverse spin polarization of the system <s> (22) 
in the state with single spin-plasmon excited as function of its 
momentum q at different values of r s . 



operators. 

To relate these quantities to the plasmon operators 
(8), we employ the unitary transformation, connecting 
operators of electron-hole excitations in noninteracting 
system a p + q iY ,a py with operators of plasmons Qq. For 
this purpose, we also need explicit expressions for 
operators of distorted single-particle excitations in the 
system with Coulomb interaction (similarly to the work 
[27]). Denoting creation operator for such single-particle 
excitation with total momentum q and energy 
£p+q,x ; — £py as Vpqyy', we can use the equation-of- 
motion method in the RPA and, similarly to the case of 
plasmons, get the explicit expression for it: 



+ _ j jYY' tt ' + 

VpqYY' ~ 2^ U pp'q a p'+q,r' a P'r' 



(24) 



prr' 



where 



tjYY I* 

pp'q 



(fpy\fp+q,Y')(. n py ~ n p+q,y') 



S 1 



(fp'+q,T' l/p'r) 



(25) 



?p'r-t, 



p+q,i 



£py + %p+q,y' + ^ 



The expressions (8), (16), (24), and (25) establish the 
unitary transformation at given q from operators of 
electron-hole excitations in noninteracting system 



a p+q,y ,a v>y to operators of excitations in Coulomb-inter- 
acting system: plasmons Qq and single-particle excita- 
tions *7pqyy'. We can easily derive the inverse 
transformation of the form: 



d p+q,y' a PY 



p'rr' 



(26) 



According to (26), we can represent the operators (11) 
and (23) of charge- and spin-density waves in the form: 



p+ = SN*n(</,CV)Q+ +/ 5+ 



st=SN*II 5 (q,CV>Qt + s; 



(27) 



(28) 



where the parts Pq and S* are the contributions of sin- 
gle-particle excitations and are dynamically independent 
on plasmons. Here, along with the usual charge-density 
susceptibility (21), the crossed spin-density susceptibility 
of the helical liquid [11] has been introduced: 



PYY' 



(29) 



x {fpy \ <7 \fp+q,y') 



co + % p+ qy + iS 



The formulas (27)-(28) show us that the average 
values of Pq and s+ in any state with a definite number 
of plasmons vanish (similarly to the mean value of 
coordinate or momentum in the simplest quantum 
harmonic oscillator). However, we can calculate 
their mean squares in the n q -plasmon 

state |n q ) = [(Q q )"V(V)~ 1/2 ]l 0 >: 



(PqPq) = <»q \PqPq\ «q> ~ (0 \PqPq\ 0), 



(30) 



(s^Y) = (n q \s^) + \ n q ) - (0 \s^) + \ 0) (31) 

(only the in-plane transverse component s ± of the spin 
s is nonzero upon averaging). Since we are interested in 
plasmons only, we have subtracted the vacuum fluctua- 
tions of these quantities in the ground state |0). 

Using (27) -(28), the mean squares of charge- and spin- 
density wave amplitudes (30)-(31) can be easily calcu- 
lated: 



(PqPq) 



nqS 2 \NqTl{q f Q q )\ 



5 q 1 (5 q 1 ) + )=n q S 2 |N q n s 1 (^^)| 2 . 



(32) 



(33) 
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We can normalize the amplitudes to obtain dimen- 
sionless quantities A p (q) = [{pqP^/riqSp] 1 / 2 and 

Ml) = [< 5 q iSqY)/nqSp] 1/2 (p=$/4n is the average 
electron density), plotted in Figure 5. As seen, the 
amplitudes of charge- and spin-density waves are close 
quantitatively at moderate momenta and any r s . 

In [11], it was shown that, due to spin-momentum 
locking, electron density and transverse component of 
spin obey an analogue of "continuity equation". It 
requires that 

Q q A p {q) = 2v ¥ qA s {q). (34) 
Our results are in agreement with this equation. 

6 Spin-plasmon scattering 

Nontrivial internal structure of plasmon in Dirac elec- 
tron gas can reveal itself in a process of its scattering on 
external potentials or impurities. For the case of spin- 
plasmon, it is also interesting to consider its scattering 
on magnetic field, acting on electron spins. Since we are 
interested in spin effects here, we will consider only in- 
plane magnetic field, affecting only spins of electrons 
and not their orbital motion. 

Hamiltonian of interaction of electrons in helical 
liquid with external electric U(r) and magnetic H(r) 
fields are, respectively, 



H e = -e ^ Up'-p(fp , Y'\fpy}al> y ,a v > y , 



pp YY' 



H m = fiB ^ Hp'_ p (/* p y \cr\fpy}a*ya P y, 



(35) 



(36) 



PP YY' 



where £/ q and H q are Fourier components of external 
electric and magnetic fields, ^ B = el2mc is the Bohr 
magneton. 

To calculate the probability of elastic scattering of the 
spin-plasmon with initial momentum q to the state with 
the final momentum q', we can use the Fermi golden 
rule, corresponding to the Born approximation. Thus, 
the differential (with respect to the scattering angle) 
probabilities of this scattering on electric and magnetic 
fields can be presented in the following form: 



dw e i i ? i . 1 9 

^ = |V q -_ q | |* e fo,0)| , 



dw m 



H, 



q-q 



*m(4,0)| 



(37) 



(38) 



where 6 is scattering angle (i.e., the angle between q 
and q'); q is absolute value of both q and q'. Here we 
have introduced electric O e (#, 6) and magnetic ® m {q, 6) 
form-factors of spin-plasmon. 



Using (35)-(36), we can calculate the form-factors 
explicitly: 



pyy'r 



(39) 



(40) 



X Cpq Cp+q'-q,q (fp+q'-q,r \fpy ) ) 
X (^pqO ( n PY ~ n p+q>')' 

*m{4,8) = imj-£r J2 ((fp + qV M/p+qr) 
Y 4 pyy'r 

xCpJq - Cp +q ,_ qq (f p+ q'_q T \cr\fpy)) 

X ( ( "pq') ( u vy ~ n p+q'Y')- 



Here ^ = dQq/dq is the derivative of spin-plasmon 
dispersion law. It is convenient to project the vector O m 
on directions, parallel and perpendicular to the initial 
plasmon momentum q to get <j)j^ and O^, respectively. 

We consider only the form-factors, revealing the spe- 
cifics of spin-plasmons, instead of the whole differential 
probabilities (37)-(38), also dependent on the form of 
external field. In Figure 6, angle dependencies of 

i 1 2 

squared modulus |O e (^ #)| 2 , <E>jL(4,0) and \<&^{q,6)\ 2 

of the form-factors at q = 0.6p F are plotted; these angu- 
lar distributions are normalized to unity. In their calcu- 
lation, we have used the single-band approximation, 
since undamped plasmons consist mainly of interband 
transitions, according to Figure 2. Results for other 
values of plasmon momentum are qualitatively the 
same. 

In the case of forward scattering with zero momentum 
transfer (at 6 = 0), the external electric field probes the 
total charge of plasmon, which is actually zero. Thus, 
the corresponding electric form-factor O e (#, 0) = 0. 
Similarly, forward scattering on magnetic field probes 
the total spin of spin-plasmon, which is directed perpen- 
dicularly to q. Therefore, ^{q, 0) = 0 and 0) ^ 0. 
As for backscattering of plasmon, it is strictly prohib- 
ited, as for individual massless Dirac electrons [28]. 

As seen in Figure 6, the form-factors demonstrate two 
side lobes, rather sharp at small r s , which can be consid- 
ered as a consequence of sharp peaking of the plasmon 
wave function (see Figure 3). At large r s , these lobes are 
much broader. 

In the process of spin-plasmon scattering on some 
configurations of electric or magnetic fields, we can 
expect interplay between spatial structure of these fields 
and that of the spin-plasmon. Angular distribution of 
the scattered plasmons will incorporate both of these 
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Figure 5 Amplitudes of charge- and spin-density waves. Amplitudes of charge- and spin-density waves. Normalized amplitudes of charge- 
density A p (a) and spin-density A s (b) waves in many-plasmon state as function of plasmons momentum q at different values of r s . 



factors, according to (37)-(38). Controlling the overlap 
of maxima of external potential and spin-plasmon form- 
factors, one can manipulate spin-plasmon scattering. 

Magnetic or nonmagnetic impurities can also create 
the external field. If the characteristic radius of the 
impurity potential is R, the transferred momentum will 
be limited by R' 1 by the order of magnitude. If qR » 1, 
the plasmon scattering on impurities will be suppressed 
due to proximity to the regime of forward scattering. In 
the opposite limit qR « 1, plasmons will be effectively 
scattered on considerable angles. 

7 Conclusions 

The properties of collective excitations (plasmons or 
spin-plasmons) in 2D gas of massless Dirac particles 



were studied. Two physical realizations of such systems 
were considered: electron gas in graphene and helical 
liquid on the surface of topological insulator. Quantum 
field-theoretical formalism for comprehensive descrip- 
tion of spin-plasmons as composite Bose particles in the 
RPA was developed. Internal structure and wave func- 
tion of spin-plasmons were studied. 

Signatures of spin-momentum locking in helical liquid 
were considered. In particular, it was shown that excita- 
tion of a spin-plasmon induces the total nonzero spin 
polarization of the system. Moreover, coupling between 
charge- and spin-density waves, accompanying a spin- 
plasmon, was demonstrated. It was shown that ampli- 
tudes of both of these waves are close by the order of 
magnitude for spin-plasmons of intermediate momenta. 
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The results of this work can be confirmed by experi- 
ments involving spin-plasmon excitation on the surface 
of topological insulator and independent measurements 
of charge and spin wave amplitudes (one of experiments 
of this type was proposed in [11]). The similar effect of 
coupling between charge and spin appears in electron 
gas with spin-orbit coupling [29], but amplitude of spin 
wave is considerable less in this case than that of charge 
wave for experimentally relevant parameters. 

Elastic scattering of spin-plasmons in helical liquid on 
electric and magnetic external fields is considered. 
Angular distribution of scattered spin-plasmons depends 
on both the shape of the potential and the form-factor 
of the spin-plasmon, revealing its internal structure. It 
was shown that, due to the form-factor, the scattering 
occurs into two side lobes, while forward and backward 
scattering is suppressed. One can use this fact to manip- 
ulate spin-plasmon scattering via interplay between plas- 
mon form-factor and shape of the external potential. It 
can be also concluded that scattering of spin-plasmons 
on long-range impurities should be very weak. 

Coupling between charge and spin waves, demon- 
strated in this article, can be used for realization of var- 
ious spintronic devices. One can perform controllable 
focusing of spin-plasmon waves and thus create regions 
with high spin polarization. Spin-polarized electrons, 
accumulating in these regions, can diffuse to adjacent 
electrodes and be used to drive spin currents (similarly 
to [12]). 

The quantum field-theoretical approach, presented in 
this article, can be used for theoretical description of 
an influence of various external factors (impurities, 
external fields) on plasmons in Dirac electron gas. The 
obtained explicit expression for plasmon wave function 
allows to derive and solve Hamiltonians of plasmons 
interacting with such external fields. In particular, the 
plasmon form-factors (39)-(40) can be used to con- 
struct matrix elements of plasmon interaction with 
external electric and magnetic fields. The problems of 
2D spin-plasmon optics, based on manipulations by 
inhomogeneities of the system and 3D environment, 
can be solved. Also the properties of hybrid modes 
(plasmon polaritons [30], plasmon-phonon modes 
[21,22], plasmon-hole modes-plasmarons [31]) can be 
studied using this approach. 

The classical electrodynamic approach based on Max- 
well equations and response functions cannot describe 
quantum effects, arising when individual plasmons are 
emitted and adsorbed. Therefore, the quantum field-the- 
oretical approach, presented in this article, should be 
especially useful for the problems of quantum plasmo- 
nics, which seems to be rather feasible in graphene- 
based structures [23]. 
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